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Abstract. The one-dimensional Ising model is easily generalized to a genuinely 
nonequilibrium system by coupling alternating spins to two thermal baths at different 
temperatures. Here, we investigate the full time dependence of this system. In 
particular, we obtain the evolution of the magnetisation, starting with arbitrary initial 
conditions. For slightly less general initial conditions, we compute the time dependence 
of all correlation functions, and so, the probability distribution. Novel properties, such 
as oscillatory decays into the steady state, are presented. Finally, we comment on the 
relationship to a reaction-diffusion model with pair annihilation and creation. 



PACS numbers: 02.50.-r, 75.10.-b, 05.50.+q, 05.70.Ln 



Complete Solution of the Kinetics in a Far- from- equilibrium Ising Chain 



2 



Introduction: With their connections to both fundamental issues of statistical 
mechanics and applications to a range of disciplines, nonequilibrium many-body systems 
have received much attention recently (see e.g. pQ and references therein). Despite these 
efforts a comprehensive theoretical approach is still lacking: As yet, there is no equivalent 
of the Gibbs ensemble theory for nonequilibrium systems. Consequently, much of 
the observed macroscopic properties of such systems are sensitive to the underlying 
microscopic dynamics, in contrast to systems in thermal equilibrium. In particular, 
most progress in this field is made by studying paradigmatic models [1, with a master 
equation governing their evolution. In this context, exact solution of simple models 
are very valuable (but also very rare), as they can be used as milestones to develop 
approximate/numerical schemes and to shed light on some general properties of related 
models. The Ising model is a good example with a venerable history QUE]- Recently, 
an interesting generalization of it was studied .3. : a kinetic Ising chain in which spins at 
alternating sites are coupled to thermal baths of two different temperatures. As a result, 
at long times, this system reaches a stationary state with inherently nonequilibrium 
properties, e.g., a non-zero heat flux through the system[3]. Subsequently, all correlation 
functions were computed exactly, so that the full stationary distribution is known [3]. 
Various other versions of two- or multiple-temperature models have been studied jSJ. 

In this letter we show that the dynamic aspects of this model are also accessible. 
As illustrations, we present the complete solution of the time-dependent magnetization 
and two-spin correlation function. These quantities can then be exploited to compute 
all other correlation functions, so that the time-dependent probability distribution is 
also (at least formally) known. The key behind the solvability of this system lies in 
the simple structure of its Glauber-like kinetics, so that the usual BBGKY hierarchy jZj 
decomposes into a closed set of linear equations for each iV-spin correlation function. 

Model specifications: We consider an Ising spin chain defined on a ring of L sites. 
For simplicity, we choose L to be even and denote a spin at site j by o~j (which 
assumes values ±1). The spins interact ferromagnetically via the usual nearest-neighbor 
Hamiltonian: TC = —JJ^jCjCj+i (J > 0; the anti-ferromagnetic case can be accessed 
by a gauge transformation). Next, we endow the system with a Glauber-like dynamics, 
but couple spins on the even and odd sites to reservoirs at temperatures T e and T , 
respectively. For T e ^ T a , this dynamics violates detailed balance [SHHE] and leads to a 
nonequilibrium stationary (but probably non-Gibbsian) state. Denoting a configuration 
(o~i,o~2, ■ ■ ■ of our system by {a}, we implement this dynamics through a master 
equation for time-dependent probability distribution P({a},t): 

dtP({a}, t) = £ [W ({a}; {a}) P({a}, t)-W ({a}; {a}) P({a}, t)] (1) 
M 

with transition rates 

^({<n;M) = £~ 

3 Z 

Here, jj is 7 e = tanh (2 J/k b T e ) for even j and 7 Q = tanh (2 J/k\,T ) for odd j. For 
convenience, the overall factor of 1/2 is chosen so that all decays follow a simple e~'-law 
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in the J = limit. 

Our goal is to compute all correlation functions (cr^ . . . <Tj n ) t = °h • • • 
. . . o~j„P({o~}, t) and to represent the complete solution for P({a}, t) by the relation 

2 L P({a},t) = l+^2cr i (a i ) t +^2a j a k (a j a k ) t + cr j a k ai(a j a k a l ) t +. . . (3) 

i j<k j<k<l 

Recently, the stationary distribution, P({a},t = oo), was found in this manner ^j. We 
will first present the solutions for (<jj) 4 and (o~jO~k)t and then, in terms of these, provide 
expressions for the other correlations. 

For Ising chains in contact with only one thermal bath, it is well known that a gauge 
transformation (changing the sign of every other spin) relates a system coupled to T < 
to one coupled to T > 0. Here, it is clear that such a transformation is applicable if the 
signs of both T e and T Q are changed. Thus, we will investigate explicitly two cases: one 
when both T's are positive and the other, when they are of opposite signs. As noted in 
jHj, quite unusual properties arise in the latter case. Here, we will find similar oscillatory 
behaviour, but in the time domain. 

The time- dependent magnetization: The single-spin function (o-j) t is, of course, just 
the t-dependent magnetisation at site j, which we denote by m,j(t). With the master 
equation (JJJ, the equation of motion of the local magnetisation reads 

J™M = ^ K-x(t) + m j+1 {t)} - m,{t). (4) 
Given any initial rrij(0), the full t-dependent magnetisation is just 
m j(fy = J2 M jk(t)m k (0) , 

k 

where Mjk(t) is the "propagator". In equilibrium (jj = 7), Glauber j2] obtained 
Mjk (t) = e~* Ik- jiffy , where J n (t) denotes usual modified Bessel function [B]. Though 
our system is not in equilibrium, we exploit this result by defining a modified 
magnetisation, nij (t) / ■Jfj, which allows to reduce (J1J) to the Glauber case with 
7 replaced by a = y/'jejo, the same parameter that enters into the steady-state 
correlations functions in [Slllj. In other words, provided T e ,T Q > 0, we can associate 
our system with an equilibrium one, coupled to a bath with temperature T eff , given by 



tanh [2J/k b T cS ] = ytanh [2J/kj,T e ] tanh [2J/k b T Q ]. To be precise, we have 

M jk (t) = e-\Pi-I k -j(act) ; a = Jfju, (5) 



which indicates that mjify suffers (linear combinations of) exponential decay similar 
to the equilibrium case. The more interesting case involves baths of opposite signs. 
Then, we can either rely on analytic continuation of 7 n (at) to pure imaginary a or 
solve equation (jlj) explicitly. The result involves oscillations with a simple exponential 
envelope: e _t . As an illustration, if the initial magnetisation is homogeneously m, then 



rrijit) 



me 



-t 



cos(|a| t) + -r-A sin (I a 1 1) 
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Interestingly, the frequency of the oscillations increases as the T's are lowered. Such 
remarkable properties can perhaps be traced to a mild form of "frustration" , arising from 
the competition between the two baths. While the effects of the positive T reservoir is 
to align spins with its neighbours, the other bath struggles to "anti-align" them. Other 
notable behaviours occur at the limits. If, say, T a — > oo (j — ► 0), then the spins at the 
odd sites decouple and m 2 j + i(t) decays purely by e~*. This allows us to integrate (jU) for 
the even sites: m 2 j(t) = m 2 j(0) e~* + r y e te~ t {m 2j _i(0) + m 2 j + i(0)} /2. At first sight, it 
may seem surprising that the effects of the neighbours linger longer. However, this aspect 
is due entirely to the details of the dynamics here (random sequential update based on 
the average spins of the neighbours). At the other extreme, there is qualitatively new 
behaviour only when both T's vanish. As expected, the uniform component of the initial 
magnetisation survives. (As a reminder, note that if the T's — > 0_, only the staggered 
component remains.) 

Equal-time correlations: Next we turn to the time-dependent two-point correlation 
function: ((Tj<Tk)t ~ (°~j)t(o~k)t- In most studies of the Ising chain, the second term is 
typically neglected, since there is generally no spontaneous magnetisation. Of course, we 
have the result for nij (t) from above and will focus only on (o~jO~k)t- The transformation 
above can still be exploited here, but it is not compatible with the "boundary condition" 
{o~kO~k)t = 1- Nevertheless, we are able to use the method of images [2] to find the 
general solution, which is rather involved and will be presented elsewhere [9 . Here, 
let us illustrate the results by restricting ourselves to a simpler case, namely, one with 
(period-2) translationally invariance (as in [3|). Then, we need to consider only four 
functions (of one variable: k — j), namely, the correlation between spins at two even 
sites, two odd sites, and one of each. We denote these by cf^(£) = (c^c^+n))*) 

c 2n(t) = (c2£-102l-l+2n)t, c|£-l(*) = {o 2l 2S ^ 2n ^ t and C° 2n _ x (t) = {o 2 i> +1 (y 2 i +2n ) t . Of 

course, (o-jO- k ) t = (akO~j)t, so that the first pair are even in n, and the last two are 
related by c e ° 2n+l (t) = <^_i(t). Thus, there is no need to study n < cases. Finally, 
we have the boundary condition (BC): c[f = Cq° = 1, the main source of complication 
here in comparison with the analysis for mj (t). 

From the master equation (JU). we find that they satisfy (for n > 0) 

^ c 2n = — 2c 2n + — [c 2n _i + c 2n+1 + c 2n _ 1 + c 2n+1 ] , (6) 

— r°° - 9r°° -\- — \r eo -\- r eo 4- r oe 4- r oe 1 (7\ 

^ c 2n — zc 2n ^ 2^ 2n ~ 1 2n+l 2n ~ 1 2n+1 ^ ' ^ ' 

^ C 2n-1 = — 2c 2n _ 1 + — [C 2n + C 2n _ 2 ] + — [C 2n + C 2n _ 2 ] , (8) 
^ C 2n-1 = — ^ C 2n-1 + ~^[ C 2n + C 2ra-2] + ~7f[ C 2n + C 2n-2] ■ (9) 

These simplify considerably, since the combinations 7 e c 2 ° — 7 c|^ and d^-i ~ c 2n-i 
decouple and just decay with e~ 2 * from their initial values. Meanwhile, the other 
combinations, 7 e c 2 ° + 7 c|^ and c 2 ^_i + c° 2n -\ are coupled, but the quantities (n > 0) 

a 2n (t) = i [ le cZ (t) + lo cZ (t)} ; a 2n _ x (t) = | (t) + c° 2 U (t)} (10) 
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allow to reduce equations (JHHHl) to a single equation: 
d_ 

dt 



^-cij = -2aj + a[aj_i + a j+1 ] , j > , (11) 



with the BC 

fl„(*)=7; l = (le + lo)/2. (12) 
Now, equations (|11I12|) are precisely those encountered by Glauber [2], the only 
differences being at, 7 instead of 7,1. An immediate consequence is the steady-state 
result, which takes the form a k (t — > 00) = 7^0 i ^0 = tanh [J/k b T e jf] , in agreement 
with those in Refs. [3J Hj (A in = Uq here). As for the complete solution with 
arbitrary initial correlations (cjCj^q, we could simply rewrite Glauber's solution here. 
Instead, let us illustrate how to derive a more compact form for the time dependence, 
in a simple example: {ujCk^jio = (i-e., uncorrected initial spins if rrij (0) = also). 
Then, 7 e C2° — J c 2n anc ^ c 2n-i — c 2n-i ( n > 0) simply remain zero for all time, so that 
c 2n°° = a 2nho,e and c^^i = «2n-i/«- To see how afc builds up to the steady-state 
value, we exploit the Laplace transforms: afc(s) = / e~ st cik(t). Condition (fT2"J) leads 
to clq(s) = 7/s , while equations (fTTj) can be solved by an Ansatz similar to the one 
in Ref.|31E|: a k (s) = A(s)u(s) k , k > 0. Inserting these into (jlEJ), we find (A ^ 0) 
+ u — (2 + s) /a. Meanwhile, ct (s) leads to A(s) = 7/s. As expected, we need only 
u = uj (s = 0) for the steady state, since the singularities of lo lie at s < —2(1 — a) 
and the pole in A controls the t — > 00 limit. For finite t, using properties of Laplace 
transforms [H] we arrive at a simple result: 

7 f 2t dr _ 

(crjcrkfr)t = ^V^li _ k \ J —e T I\ j -k\( y ar). (13) 

With this expression, we can study long-time behaviors of these correlations. Like 
the case for rrij (t), the leading decay (towards their steady-state values) is monotonic: 
£-3/2 e -2(i-a)*_ j n fa e case f or rn(t), if T e T a < 0, a turns pure imaginary, and we 
find oscillatory behavior, damped by i _3//2 e -2 * [§]. Other unusual properties emerge 
at certain limits. If, for example, T Q — > 00 (T e finite), then a,7 G — > 0. Starting 
with our initial condition, only the nearest and next-nearest neighbour correlations 
will build up to non-vanishing values |lj: c\{t) = (7 e /4) (1 — e~ 2 ') and cf(t) = 
(7 2 /8) [1 - e- 2t (1 + 2t)\ (c° 2 °{t) = 0, of course). A curious limit is T e = -T in which 
7 = 0, so that (cTjOfc^oo = 0. Thus, an initially uncorrected state never succeeds in 
building correlations. We caution that this is a somewhat singular example, as initital 
correlations are expected to survive with e~ 2t tails. Finally, we may consider the most 
extreme case: T e o — > +i _. Assuming their magnitudes are unequal in the limiting 
process, then all correlations are suppressed by 7 = O (exp {— 2J/kbT > }), where T> is 
the bigger of T e , \T Q \. 

General two-point correlation functions: Next, let us ask how a spin on site k at 
time t is correlated with the spin on site j at a later time t + t'. Following Glauber [2] 
again, we define 

c hk {t';t)= o-' j0 - k V({a'},t + t'\{a},t)P({a},t) (14) 

{<?'},{*} 
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where V({a'},t + f\{a},t), is the probability to find the system with configuration {a 1 } 
at time t + t' conditioned on the configuration being {a} at t. Being the propagator for 
the entire system, V can be represented as a sum of terms, each of which involves the 
evolution of iV-spin functions: (a^a^ . . . o~j N ) t . For our purposes here, we need only the 
first two terms: 2 L P({cr / }, t'|{er}, t) = 1 + J2k,e o~' k aeM k e (£') + ... , and arrive at 

^,k{t'] *) = E M ^ (0 to**)* • ( 15 ) 

The interpretation of this formula is clear: All correlations present at time t will be 
"propagated" by M over the time delay t' and summed accordingly. 

As an illustration, we apply (JT5j) to compute the autocorrelation function A 2 k,e{t) = 
C2k,e{t,0) — m2k{t)mt{0) for a homogeneous system with initial magnetisation m (and 
a ^ 0). The result is A 2h/ (t) = [^f e + ^f + (-l)'(V^ - e^h^at). 

This simple example shows that the amplitude of the autocorrelation function alternates 
with the parity of site I and depends on the temperatures. 

Higher correlation functions: Finally let us turn to equal-time correlations of N 
spins: (cr^ . . . o~j N ) t . For the specific case of an initially uncorrelated system with zero 
magnetisation, all functions with odd N vanish, of course. For N = 2n, we show that 
they can all be expressed in terms of two-spin functions (JT3J). 

This program is achieved by taking advantage of recent formal results obtained by 
Aliev JU]. He considered a completely general version of the kinetic Ising chain, with 
1 — £7j (cjOj-i + diCTj+i) /2 instead of [1 — 7j5j (<7j_i + o'j+i) /2] in the spin-flip rates (j2J. 
Such a model would correspond to a system with not only arbitrary nearest-neighbour 
couplings (Jfc^+i), but also a separate bath (at T k ) for each spin! In the absence of initial 
magnetisation and correlations, Aliev was able to show that the generating function for 
all correlations, ^({rj} ; t) = (1 + rjjO-j)^ , where the 77's are Grassmannian variables 

[TT] . is given by = exp (j2j<k r lj r lk (°~jO~k)t)- Using the notation of Pfaffians \H\ IT2"]. 
we can expand ^ and arrive at an expression for 2n-point function (ji < j'2 < ••• < 32n)'- 



{o h a h . . . o- hn _ x o hn ) t = 22 — -, — (o- Ul) a U2) ) t . . . (^ 7r(2n _ 1) ^ 7r(2 „ ) )t , (16) 

TT IV. 



where the summation runs over all the permutations tt of the indices 
{31J2, —,32n-i,hn}, with the constraint that j n ( 2 e-i) < jn{U) for each L Here, Sg7r 
is the signature of the permutation 71". For example, the general time-dependent 4-spin 
correlation is given by (i < j < k < £): (a i a j a k a e ) t = {eiVj)t{vhPl)t ~ {viVk)t{vj<Tl)t + 
{o~iO~e)t{o-jO-k)t- In the steady state, the last two terms cancel. More generally, such 
cancellations can be shown to persist for arbitrary n, so that we recover the result of 
Ref.@j: (a h . . .cr^Joc = (0^0^)00 . . . (o- i2n _ 1 (T j2n ) 00 . 

Before closing, we remind the readers that the results in this section hold only for an 
initially uncorrelated state with zero magnetisation. Otherwise, the iV-point functions 
will clearly be more complex than (|TE|) . Finally, due to equation (JHJ), we see that the 
full distribution P({a},t) for this specific nonequilibrium many-body problem can be 
constructed from (|T3j) and (JTKj). 
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Concluding remarks: In this letter we solved a stochastic Ising chain in which 
alternate spins are coupled to two thermal baths at different temperatures via Glauber 
spin-flip dynamics. We found analytic expressions for all correlation functions. If both 
temperatures are positive, both the steady state and the decays into it display properties 
similar to those in the ordinary Glauber-Ising model. If the temperatures are of opposite 
signs, qualitatively novel behaviours, such as oscillatory damping, emerge. Similar to 
the spatial oscillations in stationary 2-spin correlations, we believe their origins can be 
thought of as a kind of "frustration," where the two baths attempt to align/anti-align 
a spin with its neighbours. 

Finally, we note that our findings can be applied to studies of the dynamics of 
domain walls in this system As usual, the kinetic Ising model can be mapped onto 
a reaction-diffusion system (RDS), in which a "particle" (A) corresponds to a broken 
bond on the Ising lattice 0- The resulting RDS, in addition to symmetric diffusion 
(with rate 1), would be pair-annihilation (AA — > 00) with rate 1 + jj and pair-creation 
(00 — > AA) with rate 1 — "fj. Since |tj| < 1, we are satisfied that these rates are positive 
and, so, physical. Our two-temperature model is then mapped into an RDS with two 
different creation/annihilation rates on alternating sites. The implications of mapping 
our results into the RDS case are interesting, e.g., oscillatory damping of the density of 
domain- walls when 7 e 7 Q < 0. Further details will be published elsewhere ||. 
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